Fermi Gases in Slowly Rotating Traps: Superfluid vs Collisional Hydrodynamics 
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The dynamic behavior of a Fermi gas confined in a deformed trap rotating at low angular velocity 
is investigated in the framework of hydrodynamic theory. The differences exhibited by a normal gas 
in the collisional regime and a superfluid are discussed. Special emphasis is given to the collective 
oscillations excited when the deformation of the rotating trap is suddenly removed or when the 
rotation is suddenly stopped. The presence of vorticity in the normal phase is shown to give rise to 
precession and beating phenomena which are absent in the superfluid phase. 
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Recent experiments on ultracold atomic Fermi gases 
close to a Feshbach resonance (J. HJ. Q have explicitly re- 
vealed the emergence of a hydrodynamic regime which 
shows up in the anisotropic shape of the expanding gas. 
This behavior is the consequence of the anisotropy of the 
pressure gradients which are stronger in the tighter di- 
rections of the trap. It should be contrasted with the 
isotropic nature of the expansion of the non interacting 
Fermi gas, which instead reflects the initial isotropy of 
the momentum distribution. In Ref. |4( the anisotropy of 
the expansion of a Fermi gas has been suggested as a pos- 
sible signature of superfluidity. The analysis of Ref. 4] 
concerns, however, dilute gases far from Feshbach reso- 
nances where the system, in the normal phase, is colli- 
sionless at low temperature due to Pauli blocking |5j and 
consequently expands differently from a superfluid. If one 
instead works close to a Feshbach resonance the huge in- 
crease of the collisional rate will favor the achievement 
of the hydrodynamic regime in the normal phase even 
at low temperature [lj. At the same time the resonance 
effect is expected to enhance significantly the value of 
the critical temperature for superfluidity ||, providing 
promising perspectives for its experimental realization. 
Since the anisotropy of the expansion is compatible with 
either the collisional and superfluid hydrodynamic pic- 
tures, its experimental observation cannot be used as a 
test of superfluidity without further considerations. As a 
consequence, it is important to identify alternative effects 
which permit to distinguish between the two regimes. 

The purpose of the present letter is to show that the 
study of the collective oscillations of a trapped Fermi gas 
rotating at low angular velocity can provide a useful iden- 
tification of superfluidity. It is in fact well known that 
a superfluid cannot support vorticity, unless quantized 
vortices are created. In the following we will consider 
situations where vortices are absent. This can be en- 
sured by rotating the confining trap at sufficiently low 
angular velocities. Under these conditions the dynamic 
behavior of a superfluid is described by the equations of 
irrotational hydrodynamics: 



dt \2 M MJ \2 M J Mn 1 ' 

where V ex t is the external potential generated by the con- 
fining trap, while P and \i are, respectively, the pressure 
and the chemical potential of a uniform gas evaluated at 
the corresponding density. In the last identity we have 
used the T = thermodynamic relationship VP = uVfi. 
The above equations apply to dynamic phenomena of 
macroscopic type where the local density approximation 
to the equation of state is justified. They hold for both 
Bose and Fermi superfluids at zero temperature. At finite 
temperature, below T c , they should be generalized to the 
equations of two fluid hydrodynamics (see, for example, 
7J). Equations (|TJ and J2J have been systematically used 
in the last years to test the effects of superfluidity on the 
dynamic behavior of Bose-Einstein condensed gases || . 

Differently from a superfluid, a normal gas can support 
vorticity and in the collisional regime the equation for the 
velocity field takes the classical Euler form 
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where the time dependence of the pressure should be cal- 
culated taking into account the conditions of adiabatic- 
ity, following from the conservation of entropy. Equation 
differs from equation J5J) for the superfluid velocity 
because of the last term, proportional to the vorticity 
V A v. The hydrodynamic description holds pro- 
vided the collisional relaxation time r is much smaller 
than the inverse of the typical frequencies oj character- 
izing the dynamic phenomena under investigation, fixed 
by the oscillator trapping frequencies: lot <C 1. 

Let us first suppose that the initial state of the sys- 
tem does not contain any velocity flow (gas at rest in 
a static trap). In this case the equations for the ex- 
pansion as well as for the linearized collective oscilla- 
tions take the same irrotational form both in the su- 
perfluid and classical cases. In the following we will as- 
sume a trapping potential of harmonic shape: V c ^t = 
M {lj 2 x 2 + uj 2 y 2 + lj 2 z 2 )/2. If the potential is axially 
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symmetric (lo x = lo v = u>±), one can classify the solu- 
tions for the linearized oscillations in terms of the third 
component Hm of angular momentum. We will consider 
solutions where the velocity field is linear in the spatial 
coordinates. The solutions with m — ±2 and m = ±1 
are surface excitations of the form v cx V(x ± iy) 2 and 
v oc V[(x ± iy)z] respectively, with frequencies given by 



ui(m = ±2) = \/2uj±_ 
u>(m = ±1) = 



^ + ^ 2 



(4) 
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independent of the equation of state. Assuming a power 
law dependence for the isoentropic equation of state (P cx 
n 7+1 ), one easily finds analytic solutions also for the m = 
modes which are characterized by a velocity field of 
the form v cx V[a(x 2 + y 2 ) + bz 2 ). The corresponding 
collective frequencies are given by 



u 2 {m = 0) 



2( 7 + + h + 2)io 2 ± 



^4(7 + 1)2 W 4 + ( 7 + 2) 2 cj4 + 4( 7 2 - 37 - 2)w>2 
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Equation © reduces to the one derived in Ref. |l(J in 
the interacting Bose case (7 = 1), while in the case of 
the ideal gas (7 = 2/3) it reduces to the predictions of 
Refs. IllL I12L Il3l | and, for spherical trapping, to the ones 
of Refs. yjj- For elongated traps (co z -C w±) one finds 
io = ^/2(7 + 1)w_l and w = ^(3 7 + 2)/ (7 + 1) to z . 

The analysis of the collective frequencies is easily gen- 
eralized to tri- axial anisotropy (lo x ^ lo v 7^ lo z ), where 
one finds 3 solutions of the form v cx V(xy), v cx V(xz), 
v cx V(yz). These are the so called scissors modes 
|l5j relative to the three pairs of axes, with frequency 

' + to 2 , \Jio 2 x + lo 2 and \J to 2 + <-u z , respectively. The 

other three solutions have the form v cx V(ax 2 + by 2 + 
cz ) and their frequencies obey the equation 

c 6 -(2 + 7)^ 2 +^ + c 2 )c 4 + 
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Let us stress again that the above results hold both 
in the superfluid and normal hydrodynamic phases. In 
particular, if we excite the scissors mode by suddenly ro- 
tating the confining trap starting from the grou nd state 
configuration, like in the experiment of Ref. [lfj , the gas 
will oscillate with the same frequency in both regimes. 
The situation would be different in a dilute Fermi gas far 
from Feshbach resonances where the dynamic response, 
in the normal phase, is collisionless at low temperature 
and the behavior of the scissors oscillation exhibits differ- 
ent features with respect to the superfluid case 0, 0] . 
Although in a Fermi gas close to a Feshbach resonance the 
measurement of the collective frequencies (0J-0I around 



the ground state does not provide a direct indication of 
superfluidity, their experimental determination would be 
nevertheless very useful, providing an accurate and quan- 
titative proof of the achievement of the hydrodynamic 
regime. This is best tested looking at the surface exci- 
tations, whose frequencies are insensitive to the equation 
of state. 

We are now ready to explore the rotational properties 
of the system. The most natural procedure is to let a de- 
formed trap rotate at angular velocity f2 in the x-y plane. 
The angular velocity should be turned on adiabatically 
in order to ensure the conditions of stationarity. Fur- 
thermore the final angular velocity f2 should be small 
enough to avoid the formation of quantized vortices in 
the superfluid phase. This condition is not very restric- 
tive since the angular velocity needed to nucleate vortices 
by adiabatic increase of the rotation rate is very high , 
as confirmed exp erimentally in the case of Bose-Einstein 
condensates 

If the system is superfluid the stationary velocity field 
has the irrotational form v = aV(xy). The dependence 
of the coefficient a on the angular velocity fl has been 
discussed in Ref. [i(J and for small angular velocities re- 
duces to a = — efi, where e = {to 2 — lo 2 )/{lo 2 + to 2 ) is 
the deformation of the trap in the x-y plane. Actually, 
in the limit of an axisymmetric trap (e = 0) , the velocity 
field exactly vanishes revealing that the superfluid is not 
capable to rotate. If instead the system is normal, the 
stationary velocity field takes the rigid form v = fl A r, 
corresponding to VAt) = 2f2. Notice that, while the 
rigid rotation is correctly described by the hydrodynamic 
equations Q and @, the thermalization process which 
permits its achievement starting from a gas initially at 
rest, is not accounted for by Eq. (|3J), because of the ab- 
sence of viscosity terms. The time needed to achieve the 
rigid rotation of the gas (spin-up time) was calculated in 
Ref. |2ll | and in the hydrodynamic regime lo± r < 1 is 
of the order of (e 2 ^ 2 ^) -1 , where r is the relaxation time 
fixed by collisions and w 2 = (to 2 + lo 2 )/2. This suggests 
that, in order to reach steady rigid rotation in reasonable 
times, the system should not be too deeply in the hydro- 
dynamic regime and at the same time the deformation of 
the rotating trap should not be too small 29] . 

The fact that the velocity field is so different in the 
two regimes, gives rise to different predictions for the 
frequencies of the collective oscillations. For an axisym- 
metric trap the differences become particularly clear. In 
fact, while an axisymmetric superfluid cannot rotate and 
the frequencies of the m = ±2 quadrupole oscillations are 
given by Eq. (0J , in the collisional hydrodynamic case the 
degeneracy of these modes is lifted by the presence of the 
rigid rotation according to the law 



yy tni 
22]: 



io(m = ±2) = J 2 lo\ - fl 2 ± n . (8) 
Equation © is consistent with the sum rule approach 
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developed in Ref. [23J, yielding the result 



Au> = u>(m = +2) — uj(m 



-2) = A 



M (x 2 + y 2 ) 



(9) 



for the splitting of the m = ±2 quadrupole frequen- 
cies, where (l z ) is the angular momentum per parti- 
cle. In the case of a superfluid no angular momen- 
tum is carried by the system because of the irrotation- 
ality constraint, while in the collisional hydrodynamic 
regime the angular momentum is given by the rigid value 
(l z ) = MQ(x 2 + y 2 ) and hence Alj = 2 0. Measuring 
the splitting Alj then provides an efficient way to de- 
tect the vorticity of the gas. The splitting can be mea- 
sured by suddenly switching off the deformation of the 
trap. Immediately after, the system will feel the axisym- 
metric trap V e xt = M [oj"]_(x 2 + y 2 ) + to 2 z 2 ]/2, but will 
no longer be in equlibrium. Actually, in linear approx- 
imation the state of the system can be written in the 
form |0) + a + \m = +2) + a_|m = —2), where |0) is the 
new equilibrium state, while \m = ±2) are the m = ±2 
quadrupole states with excitation energies hcv± = to(m = 
±2). The coefficients a± are fixed by the initial condi- 
tions, including the quadrupole deformation and the an- 
gular velocity Q. One finds a + + a_ = (x 2 — y 2 ) / \fa 
and a+io + — cl-LO- = 2£l(x 2 — y 2 )/y/a, where a is the 
quadrupole strength [23j. The time evolution of the 
states |m = ±2) is fixed by the frequencies uj± and a 
simple calculation yields the result 

= 6 + t«m(Au,t/2) + 6-tan(q,qt) 
v ' b + - 6_tan(w Q t)tan(Awt/2) v ; 

where 9 is the angle of the principal axis of the gas in 
the x-y plane, loq = (d+ + u>-)/2 and b± — a + ± a_. 
In the superfluid case Aw = and one finds the result 
tan(2 6>) = (2 0,/loq) tan(wQt). In the collisional hydro- 
dynamic case, the angle exhibits an additional slow 
precession. This is best seen by taking stroboscopic im- 
ages at times t = 2im/L0Q, with n integer, at which 
the deformation of the atomic cloud is maximum. For 
such times tan(ujQt) — and Eq. ((TU1) yields the preces- 
sion law 8 = Atot/A. This precession is caused by the 
splitting Alj = 20 and is absent in the superfluid case. 
The numerical solution of the hydrodynamic equations 
confirms (see Fig.l) the accuracy of the prediction lfTU)l. 
based on linear approximation. In the numerical calcu- 
lation we have chosen the value 7 = 2/3 characterizing 
the equation of state of an ideal gas, including the most 
relevant case of a degenerate Fermi gas. The results are 
not however sensitive to the value of 7, consistently with 
Eq. (|TU|l . unless large values of e arc considered. The 
proposed experiment is similar to the one used in [24| to 
measure the angular momentum of quantized vortices. In 
that experiment the deformation of the gas was produced 
by suddenly switching on a laser field in an almost axy- 
symmetric Bose-Einstein condensates. This corresponds 
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FIG. 1: Time evolution of the angle 8 (measured in radi- 
ans and plotted modulus tt/2) after the sudden switch off of 
the trap deformation, (a) Superfluid hydrodynamics; (b) col- 
lisional hydrodynamics. The trap parameters are e = 0.2, 
f2 = 0.2lj_l, lo z /lu±_ — 0.1, while 7 = 2/3. The dashed line in 
(b) is the curve 6 = fit/2 (see text). 



to setting a + = a_ (&_ = 0) and Eq. (fTOf) reduces to 
6 = Auot/A. 

If instead of switching off the deformation we simply 
stop the rotation of the deformed trap, an other inter- 
esting phenomenon takes place that is worth discussing. 
In fact the gas, due to its inertia, will first continue ro- 
tating, but will soon feel the restoring force produced 
by the deformed confining potential, generating an os- 
cillation around the new equilibrium configuration. In 
the superfluid this procedure will excite the usual scis- 
sors mode with frequency \/2w^_. In the case of classical 
hydrodynamics, the gas will instead oscillate differently. 
Actually a remarkable property exhibited by equations 
and (|3J of classical hydrodynamics is that they ad- 
mit stationary solutions with non vanishing velocity flow 
also when the trap is at rest in the laboratory frame. 
These solutions have the form v = f2 A r + a'V(xy) 
and for small one finds a = eil. Under the condition 
O ^S> £ 2 uj±, the resulting oscillation can still be described 
as a linear combination of the two modes JSJ and will 
consist of the characteristic beating [^Ol 



8(t) = {CI/uiq) s'm(oj Q t) cos(Au7j/2) 



(11) 



Under the same conditions, also the intrinsic deformation 
S of the cloud exhibits an oscillatory behavior described 
by the law 5(t) — 5q = — (2ile/u>o) sm(u)Qt) sin(Aw/;/2), 
where 5q — (y 2 — x 2 ) / (x 2 + y 2 ) = e is the deformation 
of the stationary configuration. In Fig. [3 we compare 
the behavior of the angle 6 in the superfluid and nor- 
mal cases. The numerical results have been obtained by 
solving the hydrodynamic equations and in the collisional 
case exhibit the typical beating predicted by Eq. l(TT|l . 

Let us finally recall that the differences between the 
superfluid and collisional regimes discussed in this letter 
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FIG. 2: Time evolution of the angle 9 (measured in degrees) 
after the sudden stop of the rotation of the trap, (a) Super- 
fluid hydrodynamics; (b) collisional hydrodynamics. Param- 
eters as in Fig.0 



concern genuine macroscopic phenomena. The theory of 
superfluidity predicts also the occurrence of more micro- 
scopic quantum phenomena, associated with the quanti- 
zation of circulation and the appearance of quantized vor- 
tices. The description of quantized vortices requires the- 
oretical schemes beyond the hydrodynamic picture and 
has been the object of recent work in Fermi superflu- 
ids in various regimes, including the BCS [2^ and the 
unitarity limit 26]. Their observation, similarly to the 
case of Bose-Einstein condensates j^, could be again 
revealed by the splitting © of the quadrupole frequen- 
cies. In a Fermi superfluid, one predicts (l z ) = h/2 
for a single vortex line aligned along the symmetry axis 
|3l|. The realization of quantized vortices however re- 
quires different procedures with respect to the ones dis- 
cussed in the present work. In particular one should likely 
work at higher angular velocity, close to cu±/y/2 where 
the superfluid becomes unstable ag ainst the formation 
of quadrupole deformations |20l |27j. Furthermore, one 
should switch on the rotation of the trap in a non adi- 
abatic way in order to favor their nucleation. This pro- 
cedure has already proven to be successful in the experi- 
mental realization of quantized vortices in Bose-Einstein 
condensates [28). 

Stimulating discussions with L. P. Pitaevskii are ac- 
knowledged. 
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